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Abstract 

We study the spin-flavor transitions of neutrinos, v e — v^, v e — D^, etc., in the magnetic 
fields of a collapsing star. For the neutrino mass squared difference Am 2 ~ (10~ 10 — 10) eV 2 
the transitions take place in almost isotopically neutral region of the star, where the effective 
matter density is suppressed up to 3 - 4 orders of magnitude. This suppression is shown to 
increase the sensitivity of the neutrino bursts studies to the magnetic moment of neutrino, 
//, by 1.5 - 2 orders of magnitude, and for realistic magnetic field the observable effects may 
exist for /x ~ (2 — 3) • 10~ 14 /is (/j,b is the Bohr magneton). In the isotopically neutral region 
the jumps of the effective potential exist which influence the probabilities of transitions. The 
experimental signatures of the spin-flavor transitions are discussed. In particular, in the case of 
direct mass hierarchy, the spin-flip effects result in a variety of modifications of the ? e -spectrum. 
Taking this into account, we estimated the upper bounds on /iB from the SN1987A data. In 
the isotopically neutral region the effects of possible twist of the magnetic field on the way 
of neutrinos can be important, inducing distortion of the neutrino energy spectra and further 
increasing the sensitivity to uB. However, if the total rotation angle is restricted by A0 < tt, 
the absolute change of probabilities is small. 



1 Introduction 



At least 100 neutrino bursts from the gravitational collapses of stars in our Galaxy are already on 
the way to the EarthP]. A registration of even one of those by the existing detectors, Kamiokande 
M, Baksan 0, LSD 0, and especially, by the new installations: LVD || (which is already 
working), SNO fjj and Superkamiokande f§ (starting to operate in 1996), ICARUS g (which 
is at the prototype stage) will give unique and extremely rich information about features of the 
gravitational collapse, supernova phenomena and neutrinos themselves. Already SN1987A has 
given a lot. In general, one will be able to get information on the energy spectra (as a function of 
time) of f e , v e , as well as of neutrinos of the non-electron type: z/« v T , v T . Certain properties 
of these spectra do not depend on the model of the star thus opening the possibility to study the 
characteristics of neutrinos. New experiments are discussed to detect the neutrino bursts from 
collapses in other nearby galaxies [fTcfl . 

In this paper we consider a possibility to study effects of the neutrino magnetic moments 
using the neutrino bursts. In the magnetic fields the neutrinos with magnetic moments undergo 
spin-precession [11, 12, O], spin- flavor precession [14] or /and resonant spin- flavor conversion, e.g. 



v eh v^r [15, 16]. Previously, the applications of these effects to neutrinos from supernova have 
been considered in [O] (precession) and [jTJ, 17] (resonance conversion). 

The spin-flip effects are determined by the product, fiB, of the neutrino magnetic moment, //, 
and the strength of the magnetic field, B. The discovery of the effects corresponding to values of \x 
near the present upper bound, \x < (1 — 3) • 10~ 12 ^# |2l| (/j,b is the Bohr magneton) will imply rich 
physics beyond the standard model. From this point of view the studies of the neutrino bursts are 
of special interest since the magnetic fields in supernova can be as strong as (10 12 — 10 14 ) Gauss. 

It is reasonable to address two questions, keeping in mind that the magnetic field profiles of 
supernovae are essentially unknown, and vary from the star to star. What are the experimental 
signatures of the spin-flip effects? And what could be the sensitivity of the neutrino burst studies 
to the magnetic moments of neutrinos, or more precisely to /uB? 

The spin-flip probability may be of the order 1, if 

drB(r) ' ^ 

where the strength of the magnetic field is integrated along neutrino trajectory. Suggesting the 
field profile 

B-Bof-V, (2) 



r 



1 This comes about from the estimation of frequency of the collapses in our Galaxy as 1/30 year ^| 
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with Bq ~ (10 12 — 10 14 ) Gauss, k = 2 — 3, and rg ~ 10 km, where r is the distance from the center 
of star, one gets from (1), 

pB >^- (3) 

and p^(10~ 15 — 10~ 17 )pp. (For a constant field a similar bound is pB > (Arg) _1 , where Arp is 
the size of region with the magnetic field). These numbers look encouraging, being 3-5 orders of 
magnitude below the present bound. 

However, the estimation (1) corresponds to neutrino propagation in vacuum (no matter), where 
the precession takes place with maximal depth, Ap = 1. In this case the "optimal conditions" for 
the spin-flip are realized: the value of pB needed to change the spin-flip probability by a given 
amount AP is minimal p2| . The presence of dense matter strongly reduces a sensitivity to p. 
Forward neutrino scattering (refraction) results in neutrino level splitting: AH = Vsf, where 
Vsf is the difference of the effective potentials acquired by the left and right handed neutrino 
components in matter. Typically Vsf ~ Vq, where 

V = V2G F n « V2G F p/m N (4) 

is the total potential, Gf is the Fermi constant, n is the nucleon number density, p is the matter 
density and is the nucleon mass. The level splitting suppresses the depth of precession |l3|]: 



{2pBf + (AH) 2 i + [V SF /2pBf 



and to have Ap > |, one needs 

pB >^V SF - (6) 

In dense medium the restriction (6) is much stronger than (3). Indeed, for typical density of 
the neutrinosphere of the supernova: p ~ 10 12 g/cm 3 and for B ~ 10 14 Gauss, one gets from 
(6): p > 10~ 7 pp.- With increase of distance, the potential Vsf decreases as n oc r -3 , or quicker, 
and the restriction (6) relaxes faster than (3). At r>RQ it becomes even weaker than (3) (here 
Rq = 7 • 10 10 cm is the solar radius). 

The suppression of the matter density (i.e., approaching the vacuum condition) results in the 
increase of the sensitivity to p. This is realized in the case of the resonant spin-flavor conversion 



[16|. If the left and right handed neutrino components connected by the magnetic moment have 
different masses, mi, ni2, being also of different flavors, the matter effect can be compensated by 
mass splitting: 

Am 2 

AH = V SF -^ r = (7) 
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(resonance condition), where Am 2 = m\ — m\, and E is the neutrino energy. Now the expression 
for Ap becomes 

{2uBf 
{2nBY + (V SF - 

and in the resonance (7) one gets Ap = 1. However, since Vsf changes with distance, the equality 
Ap = 1 holds only in the resonance point, rp. Strong compensation of the matter effects occurs in 
some layer around rp whose size depends on gradient of the potential: Vsf = dVsF/dr. Indeed, 
from (8) one finds that A P > ±, if AH = AV S f < AVp, where 

AVp = 2uB (9) 

is the half-width of the resonance layer. The spatial size of the region with Ap > 1/2 equals 
2Ar R = 2(V SF )~ 1 AVr = 2(V S f)~ X ^B. Then the strong spin-flip effect implies that 2Arp is 
larger than half of the precession length, l p : 2Arp>l p /2 = n/2fiB. Substituting the expression for 
Arp in this inequality one gets 

_2(VB)2 

KR = ■ Zl. (10) 

tt\V S f\ 

This is the adiabaticity condition in the resonance |l^, O], and up is the resonance adiabaticity 



parameter. (Note that the adiabaticity parameter 7 used in [16] is related to that in (10) as 
Kp = 27/71"). For large densities the condition (10) is much less restrictive, than (6), allowing for a 
strong spin-flip effect for smaller values of fiB. For example, at B = 10 14 Gauss, p = 10 12 g/cm 3 , 
and Vsf /Vsf ~ r o ~ 10 km, one gets \i > 3 • lO -12 ^, instead of 10~ \ip- 

Here we will consider two new phenomena which suppress the matter effect and thus, increase 
the sensitivity of neutrino burst studies to the neutrino magnetic moment. The potential Vsf 
can be diminished for special nuclear composition of matter. We point out that this is realized 
in the external region of supernova, where the medium is almost isotopically neutral. Another 
phenomenon is the magnetic field twist - the change of the direction of the magnetic strength lines 
on the way of neutrinos. The field twist leads to an additional contribution to the level splitting 
AH which can compensate the matter effect. 

The experimental signatures of the spin-flavor transitions of supernova neutrinos have been dis- 
cussed previously in fT^| . The influence of the spin-flip on the relation between directional (induced 
by z/ e e-scattering) and isotropic (u e p - interaction) signals was studied. Without consideration of 
the dynamics of the propagation it was suggested that the probabilities of different transitions are 
either one or zero. The estimated number of events shows that four cases (1) no transitions, (2) 
only flavor transitions, (3) only spin-flavor transitions, (4) flavor plus spin-flavor transitions can 



3 



be disentangled by new installations (SNO, Super kamiokande) . We will consider a more general 
situation concentrating on the energy dependence of the transitions. 

In this paper we study the dynamics of propagation to estimate the values of fiB needed 
for different effects. We consider signatures of the spin-flip transitions related to specific density 
distribution in the isotopically neutral region and to possible presence of the magnetic field twist. 
The paper is organized as follows. In Sect. 2, we describe the effective potential in supernova. 
In Sect. 3 the level crossing schemes for different values of neutrino masses are found and the 
dynamics of neutrino propagation is considered. In Sect. 4, we study the dependence of the 
transition probabilities on the magnetic field profile and on the neutrino parameters. Sect. 5 is 
devoted to the interplay of the spin-flavor and flavor transitions. In Sect. 6 we consider some 
special effects of neutrino propagation and discuss the possible implications of the results. In 
particular, we describe the distortion of energy spectra of neutrinos and estimate the sensitivity 
of the studies of j/-bursts to [iB. Also, the upper bound on (fiB) will be obtained from SN1987A 
data. In Sect. 7, the effects of the magnetic field twist are considered. 



2 Isotopically neutral region. Effective potential. 

For the spin-flavor transitions, e.g., v e L — > v^r, v e R — > v^l, the matter effect is described by the 
potential |J, 

V SF = V2G F n(2Y e - 1) = V (2Y e - 1), (11) 

where Y e is the number of electrons per nucleon. The value of Y e depends on the nuclear compo- 
sition of the matter. In isotopically neutral medium protons = # neutrons), one has Y e = ^, 
and according to (11) V$f = 0. Consequently, the matter effect is determined by the deviation 
from the isotopical neutrality. 

For the flavor conversion the potential is proportional to the electron density: 

V F ~V2G F nY e = V Y e . (12) 

Evidently, there is no suppression of Vp in the isotopically neutral medium; Vf is suppressed in 
the central strongly neutronized region of the star. 



2.1 Isotopically neutral region 

The progenitor of the type II supernova has the "onion" structure. Below the hydrogen envelope 
the layers follow which consist mainly of the isotopically neutral nuclei: 4 He, 12 C, 16 0, 28 Si, 
32 S. Thus the region between the hydrogen envelope and the core, where elements of iron peak 
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dominate, is almost isotopically neutral. The deviation from the neutrality is stipulated by small 
abundances, £j, of the elements with small excess of neutrons, i = 22 Ne, 23 Na, 25 Mg, 56 Fe [23, p4[| . 
It can be written as 

(i_2n)«x;6.(i-^), (is) 

where Zi and A{ are the electric charge and the atomic number of nuclei a i" , correspondingly. 
Typically, one gets J2i & — 10~ 2 and (1 — ^) < 10 _1 , and therefore, according to (13): (1 — 2Y e ) < 
1(T 3 . The value of (2Y e - 1) equals 

0.6 — 0.7 hydrogen envelope 

(2Y e — 1) = ^ -(10~ 4 — 10~ 3 ) isotopically neutral region (14) 
-(0.2 — 0.4) central regions of star 

Consequently, in the isotopically neutral region the potential Vsf is suppressed by more than 3 
orders of magnitude with respect to the total potential Vq (fig.l). We calculated the potentials 
using the model of the progenitor with mass 15M & (where M & is the solar mass) |23|] . Two remarks 
are in order. For r < R@ the spatial distributions of the total density, are very similar (up to factor 
2) for stars in a wide range of masses: M = (12 — 35)M Q . The profiles may differ appreciably in the 
external low density regions which are unessential for our consideration. The collapse of the core 
and the propagation of the shock wave change both the density profile and the nuclear composition 
of medium. However, during the neutrino burst (0 - 10 seconds after the core bounce) the shock 
wave can reach the distance ~ several 1000 km at most. Therefore most part of the isotopically 
neutral region turns out to be undisturbed and one can use the profile of the progenitor. 

Strong convection effects in the inner and intermediate parts of the star, if exist, may result in 
an injection of the elements in the iron peak into outer layers, thus diminishing the degree of the 
isotopical neutrality. 



2.2 Properties of the effective potential 

According to (14) and Fig. 1, the effective potential Vsf nas the following features: It is positive 
and of the order Vq in the H-envelope; Vsf decreases quickly and changes the sign in the border 
between the hydrogen envelope and the 4 He-layer (r ~ R®)', Vsf is negative and suppressed by 
3 — 4 orders of magnitude in comparison with Vq in the isotopically neutral region; Vsf quickly 
increases in the inner edge of the isotopically neutral region (r ~ W~ 3 Rq), and becomes again of 
the order Vq in the central part of the star. In the isotopically neutral region the jumps of the 
effective potential Vsf ( n °t Vol) exist which are related to the change of the nuclear composition 
in the layers of local nuclear ignition. The jump may be as large as an order of magnitude. 
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The isotopically neutral region spreads from 10 _3 i?Q to ~ R@ (Fig. 1). Here the potential 
changes from Vsf = (10 -9 — 10~ 8 ) eV to zero: Vsf 0, when r — » (0.8 — 1)-R©. Below 10 -17 
eV, however, Vsf(t) decreases so quickly that no matter effects are induced. The detectable 
energy range of the neutrinos from the gravitational collapse is (5 — 50) MeV. Consequently, in 
the isotopically neutral region the resonance condition (7) is fulfilled for Am 2 ~ (10~ 10 - 0.1) eV 2 . 
This covers the range of Am 2 interesting from the point of view of the solar and the atmospheric 
neutrino problems. For the cosmologically interesting values, Am 2 ~ (4 — 50) eV 2 , the resonance 
is at the inner part of the isotopically neutral region, where the potential is still suppressed by 
factor 1/30 - 1/100. 

Let us stress that the suppression of the effective density takes place for the spin-flavor conver- 
sion only, i.e. when both neutrino components are active. In the case of the spin-flip into sterile 
neutrino v e —* v s , the effective potential equals V s = ^(3Y e — 1), and in the isotopically neutral 
region one has V s = The potential V s is suppressed in strongly neutronized central part of 
star, where Y e « i |]| 



3 The dynamics of the neutrino transitions 

3.1 Level crossing scheme 

We will consider the system of three massive neutrinos with transition magnetic moments and 
vacuum mixing (see [3C] for details). For definiteness we assume the direct mass hierarchy: mi <C 



77i2 -C ma and the smallness of flavor mixing. The energies of the flavor levels (the diagonal 
elements of the effective Hamiltonian) , H a , (a = v e ,i/ e ,Vr l ,i'r l ,v T ,u T ), can be written as 

H e = 0, 

H - = _y Q ( 3 y e - 1) - 0, 

A 2 (15) 
H„ = cos29 %%--V Y e , 

H p = cos 29 ^f-V (2Y e - 1) - 0. 

Here 9 is the e-/i-mixing angle in vacuum. The v T energy levels, H T , H f , can be obtained from H^, 
Hp by substituting 9 — > 9 T and Am 2 — > Am^. The angle cp(t) (cp = dc/>/dr) defines the direction 
of the magnetic field in the transverse plane. In (15) for later convenience, we have subtracted 
from all elements the energy H e to make the first diagonal element of the Hamiltonian (i.e., H e ) 
equal to zero. (This is equivalent to the renormalisation of all wave functions by the same factor 
which does not change the probabilities). The qualitative dependence of the energy levels H a on 
distance r for <p = and for different values of Am 2 /2£ is shown in Fig. 2. The resonance (level 
crossing) conditions read as H a = Hp (a / (3). 
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The peculiar behavior of the effective potential V$f i n the isotopically neutral region stipulates 
a number of features in the level crossing scheme. With increase of Am 2 one gets the following 
changes. (The levels for v e and v e are fixed, whereas v^- and ^-levels go up). 

For Am 2 ^10 _10 eV 2 (E ~ 20 MeV) the mass splitting can be neglected. Two crossings of the 
levels which correspond to the spin-flavor conversions v e — > v„ and v e — ► occur practically at 
the same point, where (2Y e — 1) ~ (r ~ O.8-R0), i.e. on the border between the H-envelope and 
4 He-layer. The crossings are stipulated by change of nuclear composition of the star (Fig. 2a). 

For Am 2 >lCT 10 eV 2 the spin-flavor resonances, v e — and v e — u^, are spatially splitted. 
With increase of Am 2 the v e — v„ resonance shifts to the center of star, whereas v e — v„ does 
to the surface. Moreover, second resonance of (P e — f M )-type appears in the outer layers. When 
Am 2 increases these same resonances approach each other, and at Am 2 ~ 3 • 10~ 8 eV 2 merge (a 
"touching" point of u e - and ^ M -levels) (Fig. 2b). 

For Am 2 ^3 • 1CT 8 eV 2 , there is only one spin-flavor resonance, v e — v^. The flavor resonance 
v e — lies closer to the surface. When Am 2 increases both these resonances shift to the center 
(Fig. 2c). 

The z/ T -level can cross v e - and zVlevels in the isotopically neutral region. Moreover, for cosmo- 
logically interesting values of Am 2 3 the ^-resonances are in inner edge of this region. They can 
also be in the central part of the star. 

3.2 Adiabaticity conditions 

In the isotopically neutral region the effective potentials for the flavor conversion and the spin-flavor 
conversion differ by 3 — 4 orders of magnitude. Therefore, the flavor and spin-flavor resonances 
are strongly separated in space, in contrast with the case of usual medium. Consequently, one can 
consider these crossings independently, reducing the task to two neutrino tasks. 

If 9 is small the flavor transition is essentially a local resonance phenomenon. It takes place at 
tr determined from (15): H^tr) = 0. Therefore, the effect depends on the local properties of the 
density distribution in resonance: Vf(vr) = VoY e (r) and Hp = Vf/Vf\ Tr - The flavor adiabaticity 
reads 

p H F Am 2 sin 2 26 

k r = • — — - • > 1. (16) 

H vr E cos26> y ' 

(This corresponds to the situation when at least half of the oscillation length, is obtained in the 
resonance layer, in the direct analogy with the definition of the adiabaticity parameter for the spin- 
flip conversion (10)). Using the potential Vf ~ Vq (fig. 1), we find that for neutrinos with E ~ 20 
MeV and Am 2 = 10~ 6 , 10~ 5 , 10~ 4 and 10~ 2 eV 2 , the resonance is situated at r = O.7i? , 0.5i? Q , 
O.17i? , O.O3i?0, correspondingly, and the condition (16) is satisfied for sin 2 29 > 0.5, 0.06, 0.02, 
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i(r 3 . 

In contrast, the spin-flip effects can be nonlocal even in case of variable density. Indeed, the 
helicity mixing angle, 9b, which determines the mixing between the left and the right components 
(in ultrarelativistic case) is 

tan20 B = 2flB A 2 - (17) 

T/ Am 2 v ' 

VSF ~ ~2E- 

(In medium with constant parameters 9 b determines the depth of precession, when the initial state 
has a definite chirality: Ap = sin 2 29b)- Since in general the field strength depends on distance, 
the angle 9b can be large enough in a wide spatial region provided fiB ~ Vsf 3> Am 2 /2E. For 
large B the spin-flip effect may not be localized in the resonance layer, and moreover, the resonance 
itself may not be local [ 30 1 . 

The rapidity of the 9b change on the way of neutrino, 9b = dflg/dr, determines the adiabaticity 
condition for spin-flip: 

/ / Am 2 

ne B <AH = J (2/i£?)2 + \V SF - — 

The adiabaticity parameter can be written as 

tt9 b ir(nBV-VfiB) 

here V = Vsf ~ Am 2 /2E. In resonance, V = or Vsf — Am 2 /2E, the expression (18) reduces to 
(10). Beyond the resonance k depends on the magnetic field change. In the extreme cases one has 



Am 2 



where kr is the resonance adiabatic parameter for a given point (10). According to (19) the 
adiabaticity condition is strongly relaxed, if Vsf 3> fJ>B above the resonance, and if Am? /2E > nB 
below the resonance. For the density and the field profiles under consideration these inequalities 
are fulfilled, so that the adiabaticity is the most crucial in the resonance. 

If the initial and final mixings are small: 9 B ~ ir/2 and 9 B ~0 and the level crossing is local, 
one can estimate the survival probability using the Landau-Zener formula: 



P « P LZ = e -~ KR = e 4 \ B *J, (20) 

where B\ = ttV/8^ 2 (see (10)). At the adiabatic condition, B = Ba, one gets Plz = 0.085. 
The probability increases quickly with decrease of the field: PlZ = 0.3, 0.54, 0.8, when B/Ba = 
0.7, 0.5, 0.3 correspondingly (see Fig. 3). An appreciable effect could be for Plz = 0.8, when the 
magnetic field is about 3 times smaller than the adiabatic value. 
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3.3 The precession and the adiabaticity bounds 

The propagation of neutrinos with magnetic moment in the magnetic field and matter is an inter- 
play of two processes: the precession and the resonance spin conversion. A relation between them 
is determined largely by density distribution. Let us define the precession bound for the product 
(fiB)p as 

VB P = ±V SF . (21) 

At (fiB) = (/j,B)p, the precession may have the depth Ap = 1/2 (see (6)); for (fiB) <C {pB)p the 
precession effects can be neglected. Let us also introduce the adiabaticity bound {i~lB)a using the 
adiabaticity condition (10): 

where Hsf = [&Vsf / &x]~ 1 Vsf is the typical scale of the potential change. For (fiB) > (^B)a the 
level crossing is adiabatic. Comparing (21) and (22) one finds that for fixed Hsf the precession 
bound increases with Vsf faster than the adiabaticity bound. Consequently, for large Vsf, the 
inequality (fJ.B)p » (fJ-B)A holds and the conversion dominates over precession. The suppression 
of Vsf makes the precession effect more profound. If the potential changes as Vsf oc r -3 , one gets 

(J*B) A oc 1 (»B) P oc 1. (23) 

Clearly, (^B)p < (ixB)a in the hydrogen envelope, and (l~iB)p > (^B)a in the center of star. The 
adiabatic conversion dominates in the inner parts, whereas the precession may dominate in the 
external layers of star (Fig. 3). In the isotopically neutral region one has (fiB)p ~ (/j,B)a, i.e., 
both processes may be essential (see Sect. 4). In fact, the condition (21) is not sufficient for a 
strong spin-precession effect. The spatial size of the region where (21) is fulfilled, Ar, should be 
large enough, so that the inequality (1) or Ar^l/fiB is satisfied. For r > 0.3R Q the latter bound 
is even stronger than (21). Instead of /jB, in further discussion we will give the estimations of the 
magnetic fields at fi = W~ 12 [ib- Rescaling to other values of fi is obvious. 

The adiabaticity and the precession bounds determined by the density distribution should be 
compared with the magnetic field profile of star. In fig. 3, we depict the profile (2) with k = 2 and 
Bo = 1.5 -10 13 Gauss. The strength of this field is practically everywhere below both the adiabatic, 
Ba, and the precession, Bp, bounds (fi = lCP 12 ^) which means that spin-flip effects are rather 
weak. The effects of the global field (2) may be strong, when Bq > 10 14 Gauss. However, for k = 3 
even extremely large central field, Bq ~ 10 17 Gauss, corresponds to very weak magnetic field in 
the isotopically neutral region. Local mechanisms of the magnetic fields generation due to possible 
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convection and differential rotation may produce the local field being much larger than the global 
field, as in the case of the Sun. 



3.4 Effects of jumps of the potential 

The electron density number n(2Y e — 1), and consequently, the effective potential V$f change 
quickly in the ignition layers, so that for reasonable values of \iB the adiabaticity is broken. This 
feature is reflected by peaks in the adiabaticity bound (fig. 3). To a good approximation one can 
consider these changes as sudden jumps of the potential, AVj = V max — V m \ n , at fixed points rj. 
The effect of jump reduces to jump-like change of the helicity mixing angle 9b- 

The jumps distort the energy dependence of probabilities. The position of the energy interval 
affected by the jump is related via the resonance condition to the values of V m { n and V^ax. The 
character of distortion depends on the magnitude of jump, AVsf, as compared with the energy 
width of the resonance layer AVr (AVr = 2fiB). (See a similar consideration for the flavor case 
in p5|]). If AVj < AVr, the jump leads to the peak of the survival probability with width, AE, 
determined by the width of the resonance layer: AE/E ~ 2fiB /VsF{ r j) |25f| - The height of the 
peak is given by the ratio: 

\AV R J \2iiB) 

For 2[iB ^> AVj the jump effect is smoothed: since the size of the resonance region is large, the 
width of peak is also large, but its height turns out to be suppressed. 

If AVj > 2AVr ~ fiB, the height of the peak is of the order 1, and the width is fixed by the 
size of jump: 

AE AVj 

E Vsf 

For AVsf ^ AVr, the jump suppresses the transition in the interval with minimal and maximal 
energies determined by V max and V^n. 

The most profound effect appears, when without jump the adiabaticity is fulfilled and the 
transition is practically complete. If the width of the resonance is small (2fiB < Vsf) an d if 
AVj ~ AVr, the jump induces a thin peak of the survival probability. Such conditions can be 
fulfilled in the inner part of the isotopically neutral region. 

Note that in the considered model of the star there are two big jumps of the potential 
(AVj /Vsf ~ 10) at r = 7 • 10~ 2 i?o and r = O.3i? (fig.l). Also small jumps (AVj/V 5F <l) 
exist. 

The elements of the dynamics presented here allow to understand all features of the evolution 
of the neutrino state. 
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4 Neutrino propagation in isotopically neutral region 



Let us consider for definiteness the spin-flavor transitions v e — and v e — v^. If Am 2 < 0.1 eV 2 
the level crossing takes place in the isotopically neutral region. For realistic magnetic fields (profile 
(2) with £ < 10 17 Gauss) the effect of 

v e — V/j, mixing can be neglected in the inner part of the 
star. Therefore, the spin-flavor evolution of the neutrino state starts from the inner part of the 
isotopically neutral region r« ~ 10' 3 R Q . At j*j the states coincide practically with pure helic- 
ity/flavor states. Flavor conversion, if efficient, can be taken into account as further independent 
transformation (sect. 5). Also the transitions which involve v T can be considered separately. 

4.1 Spin-flip effects in global magnetic field 

We use the profile (2) with k = 2 and take different values of Bq. For B$<2 ■ 10 13 Gauss the 
magnetic field strength is below the precession bound (fi = 10~ 12 /j,b) in the isotopically neutral 
region. Consequently, the depth of the precession is suppressed, and the transition has a local 
resonance character (fig. 4). The transition occurs only in the resonance region whose width 
is determined by AVr/Vsf = 2^B/Vsf- The efficiency of the transition depends on kr. The 
strongest effect takes place for the energies (over Am 2 ), E/Am 2 , which correspond to the weakest 
violation of the resonance adiabaticity (B(r) is closer to the adiabaticity bound (fig. 3, 4)). With 
the increase of Bq the E/Am 2 range of strong conversion expands (fig. 5). 

If i?o~5- 10 13 Gauss, the field profile is above both the precession and the adiabatic bounds, and 
the conversion becomes essentially nonlocal (fig. 4). For E/Am 2 > 10 3 MeV/eV 2 initial mixing 
is still very small: 6b ~ 7r/2, sin 2 2#b k, 0. Therefore, the neutrino state coincides with one of 
the eigenstates in medium. Since the adiabaticity is fulfilled, the survival probability is uniquely 
determined by the helicity mixing angle: P{y e — ► v e ) ~ cos 2 Ob- The precession effects are very 
small (fig. 4). Before the resonance the angle is determined from tan2#e ~ ^B(r) /Vsf{t). The 
mixing increases, when the neutrino enters the isotopically neutral region, and if /iB » Vsf the 
probability approaches P ~ 1/2 (fig. 4). After the resonance layer the angle diminishes quickly: 
tan20£ ~ 2/j.B /(Am 2 /2E), since the magnetic field decreases with distance, whereas the splitting 
approaches the asymptotic value Am 2 /E. Correspondingly, the probability goes to zero and one 
has almost complete spin-flip with P(v e — > v e ) ~ (fig. 4). 

For large E/Am 2 the resonance is in the external layers of the star, where the adiabaticity is 
broken. The final survival probability increases with E/Am 2 , approaching the value of probability 
before the resonance. As the result one gets P^l/2. At very large E/Am? mass splitting becomes 
unessential and the dependence of probability on E/Am 2 disappears. 

The bump in the survival probability P(E) in the energy interval E/Am 2 = (10 7 - 10 8 ) 
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MeV/eV 2 (fig. 5) is due to a density jump at r = (0.1 — O.5)i? (fig. 4 a). For very large \xB the 
effect of jump is suppressed according to the discussion in sect. 3.4. 

4.2 Spin-flip in local magnetic field 

Let us suppose that a strong magnetic field exists in some spatial region Ar between r\ and r 2 , 
so that the transitions take place in Ar only. For simplicity we consider a constant magnetic 
field B c . If fiB is below the precession bound, the conversion occurs in the resonance layer, and 
the transition probability is determined by The effect is appreciable in the energy interval 
determined via the resonance condition by the maximal and minimal values of the potential in the 
layer with magnetic field (fig. 6). When the adiabaticity is broken, the maximal effect corresponds 
to the resonance in the central part of the region with the magnetic field. In this case there is no 
averaging over the precession phase and smooth dependence of the probability on energy can be 
modulated. The energy dependence can be distorted also by jumps of the potential. 

If the field is strong enough and the adiabaticity condition is fulfilled, the average probability 
follows the change of helicity mixing angle 9 B {r) (17). The survival probability, as a function of 
E/ Am 2 is determined by the adiabatic formula: 

1 1 

P Ve ^ Ve = - (l + cos 26% cos 20Qj + - sin 26 B sin 26 f B cos $, (24) 
where <3? is the precession phase: 

$= r &rJ(2iiB) 2 + (V SF ~ Am 2 /2E) 2 . (25) 

Jr\ 

Here 6 B and 6 f B are the values of the mixing angle (17) at r\ and r 2 . For \iB 3> Vsf, the influence 
of matter is small and the survival probability is large for E/Am 2 > 1/(/j,B). If E/Am 2 » l/(/j,B) 
the precession has maximal depth (fig. 6). 

4.3 On the energy dependence of probability 

The survival probability as function of the neutrino energy has some common features for different 
magnetic field profiles. There are two characteristic energies: (E/Am 2 ) m i n and (E/Am 2 ) max . The 
former, (E/ Am 2 ) m i n , is determined by the adiabaticity violation in the inner part of the isotopically 
neutral region, or in the inner part of the region with magnetic field (in the case of the local field). 
The latter, (E / Am 2 ) mSLX , is fixed in such a way that for E > E mSLX the effect of mass splitting can 
be neglected. For the local field (E/ Am 2 ) max is determined by the potential at the outer edge of 
the region with magnetic field. In the case of global field one has (E/ Am 2 ) max w 1CT 10 eV 2 which 
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follows from the adiabaticity violation due to a fast decrease of fiB or V$f a t the border of the 
isotopically neutral region (fig. 3). 

The energies (£/Am 2 ) max and (E/ Am 2 ) m - m divide the whole energy interval into three parts. 

1) . Region of matter /vacuum suppression of the spin- flip: E/Am 2 < (E/Am 2 ) m \ n . For neutri- 
nos with such energies before the resonance layer the precession amplitude is strongly suppressed 
by matter effect. There is no appreciable conversion in the resonance region, since the adia- 
baticity is strongly broken; after the resonance the amplitude is suppressed by vacuum splitting: 
Ap ~ (fiB) / (Am 2 1 E). As a result, the probability of spin-flip is small. 

2) . Resonance region: (E / Am 2 ) max < E/Am 2 < (E / Am 2 ) m \ n . The resonance is in the iso- 
topically neutral region. If the adiabaticity is unbroken or weakly broken, the survival probability 
is small. Smooth energy dependence of P due to conversion can be modulated by the effect of 
density jumps, as well as by the precession effect, if e.g. the initial mixing angle is not small. 

3) . Precession (asymptotic) region: E/Am 2 > (E/Am 2 ) mSLX . Here the mass splitting (Am 2 /E) 
can be neglected, and therefore the spin-flip probability does not depend on energy. The spin-flip 
is due to precession and possible adiabatic (non-resonance) conversion. 

If the layers with magnetic field are in the inner part of the isotopically neutral region, the 
resonance effects can be more profound. Here the adiabaticity condition is fulfilled, even when 
the width of the resonance is small, 2/j.B/Vsf < 1- For the external layers with small density the 
adiabaticity holds for 2\iB/Vsf > 1 only, i.e. when the width of resonance is large. 

5 Spin-flip and other transitions 
5.1 Spin-flip and flavor conversion 

Let us consider the effects of the flavor conversion v e — in addition to the spin-flip. In the 
isotopically neutral region the positions of the spin-flavor resonance, r s , and the flavor resonance, 
r/, are strongly separated in space: r s < rj. If the flavor mixing as well as the magnetic field 
strength are sufficiently small, both flavor and spin-flavor transitions are local. The crossings of 
flavor and spin-flavor resonance layers are independent, and the total transition probability is the 
product of the flavor and the spin-flavor probabilities. 

Three neutrino states, (y e , v e , v^), are involved in the transitions. Correspondingly, one can 
introduce 3x3 matrix of probabilities S which relates the original, F° = (F°(i/ e ), F°(i/ e ), F ^)), 
and the final, F = (F(v e ), F(v e ), fluxes: 

F = SF°. (26) 
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Let r' be some point between the two resonances: r s <C r' <C tf. Then in the region r < r' 
the matrix S depends only on the spin-flavor transition probability, P s (pe — u a)'- S = S S (P S ). 
(Evidently, (1 — P s ) coincides with the survival probability calculated in the sect. 4). In the 
region r > r' the matrix S depends on the flavor transition probability, Pf{y e — v^) transition: 
S = Sj(Pf). The total matrix is the product S = S S (P S ) x Sf(Pf): 



( (1 - Pf) PfP s P f (l - P s ) \ 



(27) 



(1 - P 9 ) Ps 

V p f ^-Pf)P s (i-P/)(i-P.) J 

Let us consider the dependence of S on energy. If the flavor mixing is small (say sin 2 29 < 0.01), 
and the magnetic field is situated in the external part of the star, the energy regions of the strong 
flavor transition and the strong spin- flip effect do not overlap. The neutrinos of high energies flip the 
helicity, whereas low energy neutrinos undergo the flavor transition. The overlap of the resonance 
regions for the spin-flip and the flavor conversion takes place for sufficiently large mixing angles 
or/and if there is a strong enough magnetic field in the inner part of star: r < 0.3Rq (Fig. 7). 
In the overlapping energy region the probabilities P s and Pf differ from zero and the electron 
neutrinos undergo a double transition: v e — > — > v e (fig. 7). 

5.2 Effects of tau neutrino 

Inclusion of v T adds two new level crossings: the flavor, [y e — u T ), and spin-flavor, {p e — v T ) 
ones (fig. 2). Other crossings involving and v T do not give observable effects. Due to the 
assumed mass hierarchy the transitions with v T (the inner region of the star) can be considered 
independently from the transitions involving (outer part). Consequently, the total transition 
matrix is the product of the transition matrices in the inner and the outer regions. 

Three states, v e , v e , u T , are transformed in the inner part (fig. 2). For monotonously decreasing 
density, the order of resonances is the same as before: first - spin-flavor and then flavor one. If 
the v T resonances are enough separated, the transition matrix can be found from that in (27) 
by substitution Pf — > Pf and P s — > P' s , where Pf and P' s are the 2z/-probabilities of flavor and 
spin-flavor transitions in the inner region. (In fact one should consider 4x4 transition matrices in 
the (f e , iy e , i/p, z/ T )-basis in the inner and the outer parts and find their product.) The effects due 
to the z/j—resonances can be taken into account as the initial conditions, F l , for the fluxes at the 
border with the outer region. Using the matrix (27) with substitution mentioned above we get 

F(u e y = (1 - P' f )F°(u e ) + P'fP' s F\u e ) + P' f (l - P>)F°(v T ), 

F(u e y = (1 - P' s )F*(v e ) + P S 'F> T ), (28) 

F(^Y = F°(^). 
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Further on we will concentrate on the cosmologically interesting values 711,3 ~ (2 — 7) eV. This 
corresponds to Am 2 = (4 — 50) eV 2 , so that for E ~ 20 MeV the resonances are in the inner part of 
the isotopically neutral region: r ~ (1 - 2) • lO" 3 i?0. Here the total density is (0.3 - 3) • 10 8 g/cm 3 
and (l-2Y e ) ~ (2-4)-10~ 2 . According to fig. 3, the precession bound equals (10 13 -3-10 15 ) Gauss, 
and the adiabaticity bound is (0.5 — 2) • 10 11 Gauss. Thus an appreciable spin-flip transition implies 
very large magnetic field: B > 10 10 Gauss. In contrast, the flavor transitions can be efficient even 
for very small mixing angles: sin 2 29 > 10~ 6 . For Am 2 under consideration the spin-flip effects 
may be more probable below the isotopically neutral region. 



5.3 Transitions in the central region of the star 

The collapse and the shock wave propagation lead to strong time dependence of the density profile 
in the central region of the star. Moreover, the profile may be non-monotonous. With approaching 
the center of star the total matter density may first diminish from 10 8 g/cm 3 to (10 6 — 10 ) g/cm 3 
at r ~ 100 km, and then rapidly increase at the surface of the protoneutron star, r < (20 — 30) 
km. Moreover, in the central part the neutrino-neutrino scattering should be taken into account 
p§. Due to the non-monotonous behavior of Vsf, f° r Am 2 = (4 — 50) eV 2 the additional level 
crossings appear in the region r ~ (30 — 100) km: second flavor crossings v e — u T , and second 
spin-flavor crossings v e — v T . The order of resonances from the center is the following: spin-flavor 
(s), flavor (f), flavor (f), spin-flavor (s) — s, f, /, s. Note that now there is some part of the profile 
where the potential increases with distance, and consequently, the order of resonances is reversed: 
neutrino crosses first the flavor resonance and then the spin-flavor one. This may have important 
consequences for observations (see sect. 6.4). 

At distances (30 - 100) km the adiabaticity bound for the spin-flip equals Ba = (1 — 2) • 10 11 
Gauss. This is only slightly higher than that in the inner edge of the isotopically neutral region. 
However, the existence of such a field is more probable in the inner part; for the profile (2) with 
k = 2 and Bq = 10 13 Gauss the field is only 2-3 times smaller than the adiabaticity bound. Thus 
the magnetic moment fj, ~ 10~ 12 /i£ may give strong spin-flip effect here. 

If in the inner part all four transitions are efficient, there is no observable effect: the final state 
will coincide with the initial one. 



6 Implications 

In principle, future experiments will give information on the energy spectra of different neutrino 
species. Confronting these spectra with each other, one can get essentially model independent 
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information about possible neutrino transitions. It is convenient to introduce three types of original 
spectra: soft, F S (E), middle, F m (E), and hard, Fh(E), which coincide with original v e -, u e -, and 
^-spectra, respectively. The muon and tau neutrinos and antineutrinos are indistinguishable in the 
standard electroweak model at low energies, and will be detected by neutral currents. Therefore 
we will consider the total flux of the "non electron neutrinos", F{u ne ). If there is no neutrino 
transitions, then at the exit: 

F{v e ) = F s , F(U e ) = F mi F(v ne ) = AF h . (29) 

Let us find the signatures of the spin-flip transitions, as well as the sensitivity of the z^-burst studies 
to the neutrino magnetic moments. We estimate the values of B (at // = 10 _12 //b) needed for 
different effects to be observable. 

6.1 Sensitivity to magnetic moments and magnetic fields 

The suppression of the effective potential by 3 — 4 orders of magnitude diminishes the strength 
of the magnetic field (magnetic moment) needed for a strong spin-flip effect. The adiabatic and 
the precession bounds decrease by 1.5 — 2 and 3 — 4 orders of magnitude respectively. We take 
as the criterion of the sensitivity of the ^-burst studies to the spin-flip effects the magnitude of 
the magnetic field strength B s (fi = 10~ 12 /Ub) at which the probability of transition is ~ 1/2 at 
least for some neutrino energy interval. According to Fig. 2, in the most part of the isotopically 
neutral region (r < 0.3i? Q ) the adiabaticity bound is below the precession bound. Only in the 
external part of star the precession could be more preferable. Consequently, for r < 0.3-Rq the 
sensitivity limit corresponds to inequality 2fiB < Vsf- The latter means that mixing is small 
everywhere apart from the resonance layer and the transition is due to the level crossing with not 
so strong adiabaticity violation. The nonaveraged over the precession phase probability P ~ 1/2, 
can correspond to the value of the Landau-Zener probability Plz ~ 0.75. Therefore we will 
define the sensitivity limit B s (r) in a given point as the strength of the magnetic field, for which 
PLz( r ) ~ 0.75. According to the estimations in section 3.2, the sensitivity bound can be about 3 
times smaller than the adiabaticity bound: B s ~ Ba/3 (fig. 3). 

For the effect to be appreciable, the size of the region with magnetic field, Ar#, should be com- 
parable with the size of the resonance region. If Vsf oc r -3 , we get Are = (dV5^/dr) _1 2AVR ~ 
4fiBr /3Vsf, and in particular, for 2fiB/VsF < 0.3: Are < 0.2r. In fact, B s should be considered 
as the average field in the region Are. 

The transition takes place in the resonance region whose position depends via the resonance 
condition on E/Am 2 . Therefore for fixed E/Am 2 one can define the sensitivity limit B s in a 
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certain region of star r = r(V SF (E/ Am 2 )). For Am 2 = (0.3 - 1) • 10" 5 eV 2 (E ~ 20 MeV) which 



corresponds to the MSW solution to the solar neutrino problem [26, ^7| the spin- flip resonance 
lies at r ~ (1 — 3) • 10~ 2 -R©, and according to fig. 2, the sensitivity limit equals B s = (2 — 5) • 10 6 
Gauss. For the oscillation solution of the atmospheric neutrino problem (Am 2 ~ 10 -2 eV 2 ) [^] 
the resonance is at r ~ (2 — 3) • 10 -3 Rq, and B s ~ 10 9 Gauss. 

In the external part of the isotopically neutral region, r = (0.3 — 0.7) Rq, the precession bound 
is essentially below the adiabaticity bound. However, here for 2\iB ~ Vsf the precession length 
is already comparable with the distance to the center of star and the condition (3) becomes more 
important. At r ~ O.3i?0 one gets B s ~ (2 — 3) • 10 4 Gauss. This region corresponds to values 
Am 2 = (10 -9 — 10 -8 ) eV 2 which are interesting from the point of view of the resonant spin-flip in 
the convection zone of the Sun. 

In the range r ~ (0.7— 1)Rq the matter effect can be neglected and the sensitivity is determined 
by vacuum precession: B s > 2/Rq ~ (2 — 3) • 10 4 Gauss. This number can be compared with value 
of B needed to solve the solar neutrino problem. For strong v e L — > v^r conversion in convection 
zone of the sun (r ~ (0.7 — 1)-R©) one needs the magnetic field as large as B ~ 3 • 10 5 Gauss 
[ f29[ . In supernova, at the same distance from the center, the field of about 10 4 Gauss is enough. 
Conversely, if the magnetic field is ~ 3 • 10 5 Gauss, the i^-burst will be sensitive to \i ~ 3 • 10 _14 ^s- 

Note that in the isotopically neutral region one needs for appreciable v e i — v^r conversion 
1.5 — 2 orders of magnitude smaller field than for conversion into a sterile state: v e L — > v s r. 

Comparing the above results with those of sect. 5.3 we find that the sensitivity to /U can be 
even higher in the isotopically neutral region than possible sensitivity to \x in the central region of 
star. (Of course, the latter corresponds to large values of Am 2 ). 



6.2 Bounds on fiB from SN1987A 

The spin-flip — > v e leads to an appearance of the high energy tail in the iVspectrum at the 
Earth: 

F(u e ) = (1 - P s )F m + P s F h , (30) 

where P s is the transition probability. At small mixing angles the flavor conversion does not change 
this result (see fig. 3). Note that in (30) P s can be close to 1, in contrast with averaged vacuum 
oscillation effect which gives P < 1/2. 

The absence of the distortion of the i> e -energy spectrum, and in particular, the absence of the 
high energy tail, gives the bound on \xB as function of Am 2 . For a class of supernovae models 
the data from SN1987A allow to get the restriction P s < 0.35 under the assumption that Ps 
does not depend on energy Iffi] . In fact, at the border of sensitivity the suppression pit is rather 
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thin (figs. 5, 6) and the energy dependence can not be completely neglected. In this case the 
strongest observable effect takes place, when a position of the pit coincides with the high energy 
part of the z^-spectrum. The sensitivity limit obtained in sect. 6.1 gives an estimation of the 
upper bound on fiB(r). The value of the potential Vsf( t ) in the corresponding point r determines 
via the resonance condition the value of Am 2 . For example, at Am 2 ~ 1CP 8 eV 2 (resonance at 
r = 0.3R®) one gets the upper bound as B(Q.3R®) < 2 • 10 4 Gauss. Similarly, if Am 2 = 10~ 6 eV 2 : 
B(0.1R©) < 3 • 10 5 Gauss, and if Am 2 ~ 3 • 1(T 5 eV 2 : B(0MR Q ) < 10 7 Gauss etc.. In the case 
of the global magnetic field (2) with k = 2 these bounds correspond to the bound on the field at 
the surface of the protoneutron star Bq<10 13 Gauss (fi = 10 _12 //b). 

6.3 Transitions of the degenerate or massless neutrinos 

The mass differences Am 2 ^10 -10 eV 2 can be neglected for the neutrino energies E ~ 5 — 50 MeV. 
This corresponds to the asymptotic region in P-dependence on E/ Am 2 (fig. 4, 5). As we have 
noticed in sect. 2, for Am 2 ~ there are two level crossings v e — v„ and v e — Un at the same point, 
r ~ Rq . The crossings are induced by change of the nuclear composition on the border between the 
H-envelope and the 4 He-layer. However, in this region the effective density changes very quickly, 
and the adiabaticity condition implies very strong magnetic field B((0.7 — 1)Rq) > 3 • 10 5 Gauss 
which exceeds the precession bound. The effect is mainly due to the precession, and the field should 
be as strong as B ~ (2 — 5) • 10 4 Gauss. The probabilities of transitions v e «-> and v e <-> are 
equal and do not depend on energy. The final spectra of v e and v e equal: F(u e ) = (l — P s )F s +P s Fh, 
F(P e ) = (1 — P s )F m + P s Fh. In the case of complete transformation, P s ~ 1, the v e - and v e - energy 
spectra coincide with the initial f M (^)-spectrum: 

F(u e ) = F(u e ) = F h (31) 

For Am 2 > (10~ 9 — 10~ 8 ) eV 2 (which is interesting for the spin-flip of solar neutrinos |33|]) 
the effect of the mass splitting becomes important. The probabilities of v e — > Un and v e — > 9a 
transitions are different, and depend on energy. 

6.4 Distortion of neutrino energy spectra 

Using the matrix of the transition probabilities for the isotopically neutral region (27) and the 
initial conditions (28) we get for the final neutrino spectra: 

F(u e ) = (1 - P s )(l - P' s )F m + [(1 - P S )P' S + P s ]F h , 

F(u e ) = (l-P f )(l-pi)F s + [pi f pi s (l-P f ) + P f P s (l-pi s )}F m 
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+ [(1 - P f )P' f (l - P' s ) + P f (P s P' s + 1 - P s )]F h , (32) 
F{u ne ) = a s F s + a m F m + a h F h , 

where 

a s = P' f + P f (l - P' f ) 

a m = P f P' f P' s + {l- P' f )P' s + {l - P f )P s {l- P' s ) (33) 
an = 4-a s -a m . 

The important conclusions can be drawn from (32, 33) immediately. The P e -spectrum is a mixture 
of the middle (the original u e - spectrum) and hard components. It does not acquire soft component, 
except for the case of the non-monotonous change of the density in the central part of the star 
(see latter). The spectrum depends on the spin- flip probabilities P s and P' s and does not depend 
on probabilities of the flavor transitions. The changes are stipulated by the spin-flip effects only. 
Thus the distortion of the v e - spectrum, and in particular, the appearance of the hard component 
can be considered as the signature of the spin-flavor conversion. 

According to fig. 8, the spin-flip may result in a variety of distortions of the zVspectrum. In 
particular, when the transition probability P s is constant, the permutation of the v e - and v^- 
spectra can be symmetric, so that F{p e ) ~ (1 — P s )F m + P s Fh, Such an effect is realized if the 
spectra are in the asymptotic region of the suppression pit or in the region of strong (complete) 
transformation. The transition can be asymmetric, so that in P e -spectrum the suppression of the 
_F m -component is weaker than the appearance of F/j-component and vice versa. This is realized, 
when the spectra are at the edges of the suppression pit or in the modulated resonance region. 
Note that the typical energy scale of the modulations of the probability can be characterized by 
factor 2-3. Therefore the spin-flavor transition results in a smooth distortion of each component. 
The fine structure of the energy spectrum (AE ~ 1 — 2 MeV) can be due to the jumps of density 
or /and the field twist in the inner parts of the star (see sect. 7). 

Two remarks are in order. The change of the P e -spectrum is expected also if flavor mixing 
is large. In this case, however, one gets an energy independent interchange of the spectra with 
the transition probability P < 0.5. If the neutrino mass hierarchy is inverse, the properties of 
the transitions in the neutrino and antineutrino channels should be interchanged. In particular, 
zVspectrum will acquire the soft component. 

According to (33), the final f e -spectrum is, in general, an energy dependent combination of 
all three original spectra. However, F m component appears only if there are both spin-flavor and 
flavor transitions. Also the spin-flavor transition results in appearance of the F m component in 
final f M -spectrum. 
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# 




F(Ve) 


F{Vne) 


transitions 


1 


F s 


Ffn 


4F h 


no transitions 




2 


F h 


77 

1 m 


3F h + F s 


/; /'; f'f 




3 


F s 


F h 


3F h + F m 


s; s'; s's 




4 


1 m 


F h 


3F h + F s 


sf; s'f; f'sf: 


s'fs 


5 


F h 


F h 


2F h + F s + F m 


s'f; fs; s'fs; 


fs'f; fs'fs 



Table 1: The final neutrino spectra in the case when transitions are either complete or completely 
inefficient: Pi ~ or 1. In the fifth column we give a list of the transitions which result in a given 
final spectrum. Here / and /' denote the flavor transitions v e — v„ and v e — u T , s and s' denote 
the spin-flip transitions v e — and v e — v T , respectively; sf denotes a combination of the two 
complete transitions: first v e — v„ and then v e — etc. 

The modifications of the spectra are especially simple, when the resonant transitions are either 
completely efficient or completely inefficient: Pi = or 1 Jl7[ . In this case a complete permutation 
of the original spectra occur. Moreover, as follows from (33, 34), there are only four possible types 
of final spectra (Table I). 

According to the Table I, the flavor transitions (# 2) result in a hard ^ e -spectrum, whereas 
z/ e -spectrum is unchanged. Correspondingly, the f ne -flux acquires a soft component. The neu- 
tronization peak consists of i^ ne -neutrinos. The spin-flavor conversion only (# 3) leads to a hard 
z/ e -spectrum and to unchanged f e -spectrum. The flavor conversion in the inner part of star could 
be efficient during the early stage of the burst, so that the spectra #3 for cooling stage can be 
accompanied by the i/ ne -neutronization peak. The spin-flavor and the subsequent flavor transitions 
of v e in the inner or outer parts of the star result in the f e -spectrum coinciding with the initial 
z/ e -spectrum, F m (# 4). In this case v e has hard spectrum Fh- The same final spectra (# 4) 
appear also if in addition some other transitions take place which do not influence the neutrino 
flux originally produced as v e . The flavor transition of v e and the spin-flavor transition of v e result 
in the same hard spectra for v e and v e (#5). In principle, future experiments will be able to 
distinguish these possibilities. 

If the transitions are incomplete, the final spectra are certain energy dependent combinations 
of the above five spectra. 

The spectra can show a strong time dependence. The original spectra themselves depend on 
time: the temperatures decrease, and moreover, they decrease differently for different neutrino 
species. Also the transition probabilities change with time due to the variations of the effective 
potential profile in the inner part. 
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Interesting effect can be related to the existence of a region in which the density increases with 
distance. As we have noticed, here the order of resonances is changed: neutrinos first cross the 
flavor resonance and then the spin-flavor resonance. In this case v e can be transformed into P e ; 
the final i/ e -spectrum will contain the soft component, and moreover the neutronization peak will 
consist of V e . Such a modification of spectra can be realized, if two inner resonances s and /, are 
inefficient. They are at the surface of the protoneutron star, where the density is changed very 
quickly, and therefore the adiabaticity could be strongly broken. 



7 Effects of magnetic field twist 
7.1 Field twist. Scale of twist 

If the direction of the magnetic strength lines changes with distance, the propagating neutrinos feel 
the rotation (twist) of the magnetic field. The field twist leads to the splitting of the levels with 
different helicities by the value (see (15)) The effect can be considered as the modification 
of the effective potential: 

V SF -^V 4> = V SF + 4>- (34) 

Thus the field twist changes the level crossing picture (the positions of resonances). However, since 
the flavor and spin-flavor resonances are strongly separated due to the isotopical neutrality of the 



medium, the shift (34) hardly induces the spatial permutation of resonances [30]. 
The field twist can be characterized by the scale of twist, r^: 

IT 

r = T ' ( 35 ) 

9 

so that on the way, r^, the total rotation angle (in case of uniform rotation) equals A<p = tt. 
It is natural to suggest that the total rotation angle is restricted by 

A<p < vr. (36) 

For example, the twist appears, when the neutrinos cross the toroidal magnetic field with strength 
lines winding around the torus. In this case the maximal rotation angle is tt, i.e. the bound (36) 
is fulfilled. Let Ars be the size of the region with the magnetic field, then A(f> ~ 4>ArB, and the 
condition (36) implies 

r^>Ar B . (37) 

That is the scale of the rotation is comparable or larger than the region with the magnetic field. 
For the global field one has Are ~ r, and therefore 

r> > (0.1 - l)r. (38) 
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As we will see the bound on the total rotation angle (36) restricts strongly the effects of the field 
twist. 



7.2 Effect of density suppression 

Let us consider possible increase of sensitivity to (nB) due to the field twist. According to (34) at 
<ft = (fic, where 

k = -V SF (r), (39) 

the matter effect is completely compensated at the point r. Let us define the critical rotation scale, 
r£, at which the condition (39) is fulfilled: 

ri-~- (40) 

(Note that r£ coincides up to the factor 2 with the refraction length). At ~ the effect of 
field twist can be essential, and for < r£, it even dominates over the density effect. Comparing 
r£ with the distance from the center, r, we find (fig. 9) that r^/r > 1 for r > 5 • lO~ 2 i? , and 
<C r for r < 5 • lO~ 3 i? . In the inner parts of the star the critical scale is much smaller than the 
distance from center, e.g., at r = 10~ 3 R & one gets ~ 20 cm. The condition (38) is fulfilled for 
r ~ 5 • lO~ 3 i? . Let us stress that the isotopical neutrality essentially enlarges the region, where 
(38) is satisfied and therefore the field twist can be important. 

Suppose the equality (39) is fulfilled in some layer of size Ar (evidently, the field twist should 
be nonuniform). In Ar the spin-flip effect has a character of precession with maximal depth, and 
the sensitivity to [iB is maximal. Let us estimate \xB taking into account the restriction (36). The 
total rotation angle Acp equals 

A(f>= [ <p(r)dr = - [ V SF (r)dr, (41) 

J Ar J Ar 

and if Ar <C r: 

A(p ~ -V S f Ar. (42) 

The spin-flip probability is of the order one, when Ar ~ ir(2fiB)^ 1 . Substituting this Ar in (42), 
we get the relation: 

" B *W ' (43> 

Finally, (43) and the bound on the total angle of the twist (36) give fiB > Vsf/% which precisely 
coincides with the precession bound (6). Thus, the field twist can relax the precession bound by 
factor 1/2 at most, and to the further increase of the sensitivity one should admit the field rotation 
on the total angle which exceeds tt. 
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In the external layers (r > O.I-R0), Ar is larger than r, and the compensation (39) implies fine 
tuned profiles of <p(r) and V(r) in a wide spatial region which seems rather unnatural. On the 
other hand, if the compensation takes place in the inner part (r < 0.1R & ) then Ar <C r. However, 
here the precession bound is more stringent than the adiabaticity bound, and the compensation 
does not allow to gain in the diminishing of the field strength. In the inner part the effect of the 
field twist can be due to influence on the adiabaticity. 

7.3 Influence of field twist on the adiabaticity 

Nonuniform field rotation ((f) / 0) modifies the adiabaticity condition. The adiabaticity parameter 
equals: 

H = ^t. (44) 
\v SF + 4>\ 1 ' 

At 

4> ^ -V SF (45) 

K<f) — > 00, i.e. there is a strong improvement of the adiabaticity. The field twist results in flattening 
of the potential, so that dV^/dr ~ or — const., in some region Ar. Let us estimate the 
minimal value of the total rotation angle in this region. For this we suggest that <p = at the one 
of the edges of Ar. Then from (45), it follows that A</> ~ Vs^(Ar) 2 /2. Inserting Ar ~ ir(2fiB)~ 1 
(the condition for the strong precession effect) in the last expression, we get 

A0« , (46) 

V 2(2/i£?) 2 kr V ; 

i.e. the total rotation angle is the inverse value of the adiabaticity parameter without field twist. 
Therefore, if the adiabaticity is strongly broken, kr <C 1, one needs Acj) » ir to get an appreciable 
spin- flip effect. For restricted values of total rotation angle (36) only weakly broken adiabaticity 
can be restored. The relation (46) means that if the transition without field twist is weak (k <S 1), 
then field twist will induce a weak effect for A0 < tt. On the other hand, if the transition without 
twist is strong (kr ~ 1), then the field twist can make it even stronger but the absolute change of 
probability turns out to be always small (fig. II). 

The effect of the adiabaticity restoration due to field twist can be considered as the effect of 
precession in the region Ar with flat potential (V^ = const.). For constant B (which we suggest 
for simplicity) the transition probability equals 



p (W sin 2.L.. m2 , U Am^Ar 



(2/i£?) 2 + (V^ - Ar^£^2 



2E 



+ (47) 
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One can find from (47) that field twist leads to the distortion of the energy dependence of the 
probability in some energy interval AE located at E/Am 2 ~ (2V^,) _1 . Let us estimate the size of 
AE. According to (47) AE is determined by the resonance width AE/E ~ 4/xB/Vsf- In turn, 



2/iB can be estimated from (46) and the condition A(f> < it: 2fiB>y (ttVsf)/^ ■ So that finally, 
we get 



M>V 2 ^ = . ^ (48) 
E Vsf V HsfVsf 

If Vsf oc r ~ 3 an d Hsf oc r ) the relation (48) gives ^ oc r. The closer the region with field twist 
to the center of the star the thinner the peak. For r ~ W~ 2 Rq, one gets ^^>0.3. In external 
region of star (r > 0.1i?o), the field twist results in smooth change of the probability in a wide 
energy region, and it is impossible to distinguish it from other effects (fig. 11). In any case the 
absolute value of AP is small, unless A(p ir. 

As follows from (44) the field twist can also destroy the adiabaticity, when \ 4>\ 3> \ V\. Randomly 
twisting field inhibit the conversion in such a way. This can result in bumps of the survival 
probability that may be similar to those due to density jumps. 



8 Conclusion 

1. For neutrino mass squared differences Am 2 < 10 eV 2 (which are interesting for the cosmology, 
as well as for the physics of solar and atmospheric neutrinos) the resonant spin-flavor transitions 
{y e — ► Vp, etc.) take place in almost isotopically neutral region of collapsing star. In this region 
which extends from ~ 10~ 3 R Q to ~ R Q the deviation from the isotopical neutrality, 2Y e — 1, can 
be as small as 10~ 3 — 10~ 4 . Correspondingly, the matter potential for the spin-flavor transitions, 
being proportional to (2Y e — 1), turns out to be suppressed by 3 - 4 orders of magnitude. Moreover, 
the potential changes the sign at the inner edge of the hydrogen envelope. 

2. The suppression of the effective potential in the isotopically neutral region diminishes the 
values of (/jB) needed to induce an appreciable spin-flip effects by 1.5 - 2 orders of magnitude. 
Thus, the sensitivity of the neutrino burst studies to the transition magnetic moments of neutrinos 
increases, being of the order 10 _13 //-B for Am 2 = (10~ 8 — 1CT 1 ) eV 2 and for a reasonable strength 
of the magnetic field. In particular, for Am 2 = (10~ 8 — 10~ 9 ) eV 2 the desired values of (fiB) 
turn out to be 1.5 — 2 orders of magnitude smaller than those for a strong conversion of the solar 
neutrinos. 

3. In the isotopically neutral region the potential changes very quickly in the layers with 
local ignition (jumps of the potential). The jumps result in distortion of the energy dependence of 
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probabilities. In particular, one can expect the appearance of thin peaks in the survival probability 
if the jump is situated in the inner part of the star. 

4. Depending on the values of the neutrino parameters as well as on the magnetic field profile 
one expects a variety of the modifications of the neutrino spectra. In the case of a direct mass 
hierarchy and a small flavor mixing the main signature of the spin-flip effect is the distortion of the 
z? e -energy spectrum, and especially the appearance of the high energy tail. In general, the final v e - 
spectrum is the energy dependent combination of the original t> e -spectrum and the hard spectrum 
of the non-electron neutrinos. Another important signature of the spin - flip can be obtained from 
the comparison of the spectra of different neutrino species. In particular, the v e - and z^-spectra 
can be completely permuted. The combination of the spin-flip effect with other (flavor) transitions 
may result in rather peculiar final spectra. For example, v e may have the spectrum of the original 
v e , whereas v e may have the original ^-spectrum. The electron neutrino and antineutrino spectra 
can be the same and coincide with the hard spectrum of the original muon neutrinos, etc. 

5. The resonant spin-flip effect for the massless or the degenerate (Am 2 <lCT 10 eV 2 ) neutrinos 
can be induced by the change of the nuclear composition. Such a transition takes place both for 
neutrinos and antineutrinos in the same spatial region (near to the bottom of the H-envelope), 
resulting again in hard and equal u e - and P e -spectra. 

6. The absence of the considered effects allows one to get the bound on (fiB(r)) as a function 
of Am 2 . We estimate such a bound for SN 1987A. 

7. In the isotopically neutral region, where the matter potential is strongly suppressed, the 
effects of (even) small field twist may be important. Field twist can further suppress the potential, 
thus increasing the sensitivity to the magnetic moment. However, if the total rotation angle is 
restricted (< ir), a possible diminishing of \xB can be by factor 1/2 at most. The twist of the field 
may induce a distortion of the neutrino energy spectra. In particular, improving adiabaticity, in 
the inner part of the isotopically neutral region the field twist may lead to peaks of the transition 
probabilities. 

Acknowledgements 

The authors are grateful to E. Kh. Akhmedov for discussions and for numerous remarks concerning 
this paper. H. A. would like to thank Prof. A. Salam, the International Atomic Energy Agency 
and UNESCO for hospitality at the International Centre for Theoretical Physics. 



25 



References 

[1] S. Van den Bergh, Phys. Rep. 204, 385 (1991). 
[2] K. Hirata et al., Phys. Rev. Lett. 58, 1490 (1987). 
[3] R. M. Bionta et al., Phys. Rev. Lett. 58, 1494 (1987). 
[4] E. N. Alexeyev et al., Phys. Lett. B 205, 209 (1988). 
[5] M. Aglietta et al., Europhys. Lett. 3, 1315 (1987). 

[6] C. Bari et al., Nucl. Inst. Meth. A 277, 11 (1989); M. Aglietta et al., Nuov. Cim. A, 105 1793 
(1992). 

[7] Scientific and Technical Description of Mark II SNO Detector, ed. by E. W. Beier and D. 
Sinclair, Report SNO-89-15 (1989). 

[8] K. Nakamura, ICRR-Report-309-94-4. 

[9] P. Benetti et al., Nucl. Inst. & Meth. A315, 223 (1992), ibidem A332, 395 (1993). 

[10] D. B. Cline, preprint UCLA-CAA 107-3/94. 

[11] A. Cisneros, Astrophys. Space Sci. 10, 87 (1970). 

[12] K. Fujikawa and R. E. Shrock, Phys. Rev. Lett. 45, 963 (1980). 

[13] M. B. Voloshin, M. I. Vysotsky, and L. B. Okun, Zh. Eksp. Teor. Fiz. 91, 754 (1986) [Sov. 
Phys. JETP 64, 446 (1986)]. 

[14] J. Schechter and J. W. F. Valle, Phys. Rev. D 24, 1883 (1981); ibidem 25, 283 (E) (1982). 

[15] C.-S. Lim and W. J. Marciano, Phys. Rev. D 37, 1368 (1988). 

[16] E. Kh. Akhmedov, Sov. J. Nucl. Phys. 48, 382 (1988); Phys. Lett. B 213, 64 (1988). 

[17] E. Kh. Akhmedov and Z. G. Berezhiani, Nucl. Phys. B 373, 479 (1992). 

[18] M. B. Voloshin, Phys. Lett. B 20, 360 (1988). 

[19] J. T. Peltoniemi, Astron. & Astrophys., 254 121 (1992). 

[20] E. Kh. Akhmedov, Sov. Phys. JETP 68, 676 (1989). 



26 



[21] M. Fukugita and S. Yazaki Phys. Rev. D 36, 3817 (1987); G. G. Raffelt, Astrophys. J. 365, 
559 (1990); V. Castellani and S. Degl'Innocenti, Astrophys. J. 402, 574 (1993). 

[22] M. Moretti, Phys. Lett. B 293, 378 (1992). 

[23] S. E. Woosley and T. A. Weaver, Annu. Rev. Astron. Astrophys. 24, 205 (1986). 
[24] S. E. Woosley, N. Langer and T. A. Weaver, Ap. J. 411, 823 (1993). 
[25] P. I. Krastev and A. Yu. Smirnov, Mod. Phys. Lett. A 6, 1001 (1991). 
[26] S. P. Mikheyev and A. Yu. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985). 
[27] P. I. Krastev and A. Yu. Smirnov, Phys. Lett. B 338 282 (1994). 

[28] R. Becker-Szendy et al., Phys. Rev. D 46, 3720 (1992); Y. Fukuda et al., Phys. Lett. B 335, 
237 (1994). 

[29] see for latest discussion E. Kh. Akhmedov, A. Lanza and S. T. Petcov, SISSA Report 
169/94/A-EP. 

[30] A. Yu. Smirnov, Phys. Lett. B 260, 161 (1991); E. Kh. Akhmedov, S. T. Petcov, and A. Yu. 
Smirnov, Phys. Rev. D 48, 2167 (1993). 

[31] A. Burrows, D. Klein and R. Gandhi, Phys. Rev. D 45, 3361 (1992). 

[32] A. Yu. Smirnov, D. N. Spergel, and J. N. Bachall, Phys. Rev. D 49, 1389 (1994). 

[33] E. Kh. Akhmedov, Phys. Lett. B 213, 64 (1988). 

[34] L. B. Okun, Sov. J. Nucl. Phys., 48 967 (1990); J. Pantaleone, Phys. Lett. B 287, 128 (1992); 
S. Samuel, Phys. Rev. D 48, 1462 (1993). 



27 



Figure Captions 

Figure 1. The dependence of the potentials for the spin-flavor conversion, u e — ► (dotted and 
solid lines) and the flavor conversions, v e — > z/^, (dashed line) on the distance from the center of 
the star (M = 15M ). At r ~ O.8-R0, the spin-flavor potential changes the sign, and for r > 0.8Rq 
we depict — Vsf which coincides with Vsf f° r v e — * channel. 

Figure 2. The qualitative dependence of the energy levels on the distance from the center of the 
star, r, for different values of Am 2 /2E. The energies are defined in such a way that H{y e ){r) = 
(solid line). (The energy scale is relative). Dotted lines show a behaviour of the eigenvalues of the 
Hamiltonian near the crossing points. Also marked are the positions of the inner and outer edges 
of the isotopically neutral region, a). E/Am 2 « 0; b). E/Am 2 <3 ■ 1CT 8 eV 2 ; c). E/Am 2 >3 ■ 1(T 8 
eV 2 ; 

Figure 3. The dependence of the adiabaticity Ba (solid line) and the precession Bp (dotted 
line) bounds for fi v = 10~ 12 /iB on the distance from the center of star. Also shown is the magnetic 
field profile (2) with B = 1.5 • 10 13 Gauss and k = 2 (dashed line). 

Figure 4. Spatial picture of the {p e — > f M )-spin-flavour transition. The dependence of the 
zVsurvival probability on the distance from the center of the star for different strengths of the 
magnetic field. The field profile (2) was used with k = 2 and Bq (in Gauss): 5- 10 12 (dashed-dotted 
line), 1.5-10 13 (dashed line), 5-10 13 (dotted line), 1.5-10 14 (solid line), a). E/Am 2 = 10 8 MeV/eV 2 
b). E/Am 2 = 10 6 MeV/eV 2 . 

Figure 5. The i/ e -survival probability for the (i> e — > f M )-spin-flip transition in the global mag- 
netic field as a function of E/Am 2 for different strengths of the magnetic field. The magnetic field 
profile (2) was used with k = 2 and Bq (in Gauss): 1.5 • 10 12 (bold solid line), 5 • 10 12 (dashed line), 
1.5 • 10 13 (dashed-dotted line), 5 • 10 13 (dotted line), 1.5 • 10 14 (solid line). 

Figure 6. The i/ e -survival probability for the (p e — ► f M ) - spin-flavor transition, in the local 
magnetic field as functions of E/Am 2 for different values of the magnetic field. A constant magnetic 
field B c in the internal (1 — 3) • W~ 2 R Q was used. B c (in Gauss): 3 • 10 7 (solid line), 6 • 10 7 (dashed 
line), 2 • 10 s (bold solid line), 3 • 10 8 (dashed-dotted line). 

Figure 7. The probabilities of different transitions for the neutrino system with a magnetic 
moment and a vacuum mixing as functions of E/Am 2 . The initial state is (a) i> e , and (b) u^. The 
lines show the probabilities to find in the final state u e (solid), u e (dotted), and (dashed line). A 
vacuum mixing angle sin 2 20 = 1CT 2 and a constant magnetic field B c = 10 4 Gauss in the interval 
R = (0.1 - 1.3)i? were used. 
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Figure 8. The distortion of the z/ e -spectrum. The dependence of the product (flux) x (energy 
squared) on energy for different sets of neutrino parameters and different configurations of the 
magnetic fields. The original P e -spectrum is shown by the bold solid line. The case of the complete 
— v a transformation {F = F^) is shown by the solid line. 

Figure 9. The critical scale of the field twist over the distance from the center of the star as a 
function of r. 

Figure 10. The effects of the non-uniform field twist. The ^-survival probability of the spin- 
flip transition, v e —> z/„ , as a function of E/Am 2 for different configurations (scales) of the field 
twist. The field profile (2) with k = 2 and Bq = 5 ■ 10 12 Gauss was used. The solid line shows 
the probability without field twist. The field twist is located in the region tq — {vq + Sr), and the 
profile of twist is described by <j>(r) = (2ir/6r)(l — (r — ro)/6r); the total rotation angle equals ir. 
The curves correspond to r = 5- KT 3 i? and Sr/R Q = 1(T 4 (dotted), 5- 1(T 4 (dashed), 2.5- 1(T 3 
(dashed-dotted). 
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